function ejer1b()
    a = 0;
    b = 10;
    N = 20;
    
    x = linspace(a, b, N);

	% Coeficientes de polinomio de grado 1: f(x) = (3/4)*x - 1/2
    coef = [3/4, 1/2]; 

	% Evalua el polinomio, con dispersion normal en el eje y
    y = polyval(coef, x) + randn(1,20);

	% Devuelve coeficientes del ajuste polinomial de grado 1.
    coef_ajustados = polyfit(x,y,1);

	% Evalua el polinomio que ajusta.
    y_ajustados = polyval(coef_ajustados, x);

    plot(x,y,'r');
    hold on
    plot(x,y_ajustados,'b')
    	
endfunction
	